Characteristics of the Wave Function of Coupled Oscillators in Semiquantum Chaos 
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Using the method of adiabatic invariants and the Born-Oppenheimer approximation, we have 
successfully got the excited-state wave functions for a pair of coupled oscillators in the so-called 
semiquantum chaos. Some interesting characteristics in the Fourier spectra of the wave functions 
and its Correlation Functions in the regular and chaos states have been found, which offers a new 
way to distinguish the regular and chaotic states in quantum system. 
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I. 



INTRODUCTION 



Even though plenty of research work has been done on 
chaos in quantum system, the surprising fact is that its 
proper definition has not been clearly got now. 

From "the principle of quantum corresponding" , there 
must be counterpart of classical limit for a quantum dy- 
namical system. Percival [1] showed that the semiclassi- 
cal approximation provides a way to divide the quantum 
energy spectrum of a N-dimension conservative system 
into two parts, regular and irregular, which relate, re- 
spectively, to periodic and chaotic motions in their clas- 
sical phase space. Bohigas et. al. [2] and Berry et. al. [3] 
then developed the approach of energy-level statistics to 
describe the properties of regular and irregular quantum 
energy spectra. Meanwhile, the studies of eigenfunctions 
made by Shapiro [4] and Heller [5] give us more informa- 
tion about the regular and irregular states of quantum 
systems. 

In quantum mechanics, nonexistence of the concept of 
path or trajectory prevents us from using the definition of 
classical chaos, which is characterized by the sensitivity 
of the system to the initial conditions. There are many 
theories to study properties of the irregular behavior, or 
the so-called quantum chaos, but most researchers have 
focused on the concept of integrability in quantum me- 
chanics and hope to get an exact definition of the quan- 
tum chaos. 

In this paper, we will consider a system in which 
a classical oscillator interacts with a purely quantum- 
mechanical oscillator. This system was first considered 
by Cooper et.al. [6], which can be described by a clas- 
sical effective Hamiltonian, the expectation value of the 
quantum Hamiltonian. Faccioli et.al. in their comment 
[7] gave a more proper treatment for this problem. But 
Cooper et.al. and Faccioli et.al. had only paid their 
attention to the average value of the time-dependent oc- 
cupation number in the chaotic state. As we all know, in 
quantum mechanics, the wave function can give not only 
the static but also the dynamical information of the sys- 
tem. Recently, Hui et. al. in their paper [8] studied the 
ground-state wave function of the coupled oscillators, and 
gave out a method to distinguish the chaotic evolution in 



this system. But they did not study the behaviors in its 
excited state. Then there come some questions: Will this 
method still work in the excited-state? Could we get also 
the excited-state wave functions for the coupled oscilla- 
tors? What are the properties of the excited-state wave 
functions? 

In this paper, we will try to answer these questions, and 
then to find whether there is another way to reach the 
goal. In section II, we will use the method of adiabatic 
invariant [9] and the Born-Oppenheimer Approximation 
[10-13] to obtain the wave functions in any excited-states. 
In section III, we will use these wave functions to discuss 
the time evolution characteristics of the system, and cal- 
culate their correlation functions. Finally, we will end the 
article with some conclusions and discussions in section 
IV. 



II. EXCITED-STATE WAVE FUNCTIONS OF 
THE COUPLED HARMONIC OSCILLATORS 

In this paper, we mainly deal with a coupled system, in 
which a classical oscillator interacts with a purely quan- 
tum mechanical oscillator, described by a classical La- 
grangian: 
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where the coordinates x and A describe, respectively, the 
motion of the quantum oscillator and the classical one 
in the system. Using the Born-Oppenheimer approxima- 
tion, we can decouple the system into a classical part and 
a quantum part, each of which can be handled easily. 

From Lagrangian (1), a Schrodinger equation for a 
state with total energy i?„can be obtained: 
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-0 {A) Xn {x, A) and following the treatment in Ref. [7,i 
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i.e., using the semiclassical approximation, we obtain the 
following coupled equations for the tp (A) and Xn i^i ^) 



On defining the nth excited state of a by 
(a)"^^ \n) = 0, we can get in the coordinate rep- 
resentation, 



and 
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0„(x,G(i)) = 0, (12) 
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where 0„ [x, G (t)) is just the nth excited state wave func- 
tion, i.e., \n). Its normalized solution can be easily ob- 
tained as follows: 



where 



Xn {x, A) = Xsn {x, A) exp 
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and for an operator O, its average value in nth excited 
state is defined as 



O) = 



I dxX*nXn 



Now, it is easy to solve the Schrodinger equation, Eq. 
(jll), by using the adiabatic invariant method mentioned 
in Ref. [1,7,8], in which an adiabatic invariant I{t) and the 
time-dependent canonical lowering and raising operators 
d and d^ are introduced. In order to compare with the 
paper of Cooper et.ai, we let G = {x'^)^. Thus, d has 
the following form: 
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where H„ (x) is the nth Hermit function, defined 

as: Hn (x) = (-1)" (2"n!V^)"5 exp (x^) ^ exp {-x^). 
i^nd we take h — I and m — 1 for simplicity. 
From Eq. (|13p. we have 
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, which is just the ground 

state wave functions given in Ref. [8]. 

And it is easy to get the total energy and the quantum 
part energy: 



En = \a^ + (2n + 1) 
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Further, when we have given the fixed total energy En, 
the A and G satisfy the following equations: 
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Especially, when n = 0, Eo = \A^ 
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\A^ + Eqq. 

In order to compare with the results in Ref. [8], we 
will only fix the total system energy (then, the initial 
condition and the phase space are fixed.) in the following 
discussion. 



III. DISCUSSIONS ABOUT THE 
EXCITED-STATE WAVE FUNCTION 



Utilizing Eq. p^ . we can obtain the Fourier spectra 
(FS) of probability density operator in the excited-state 
wave functions in chaotic or regular state at the fixed 
position X. And from them, we can see whether the phe- 
nomena happened in ground-state case would be possible 
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^ plot of the suifaoe ot seotion for energy=2 7,e=l .0,n=1 and A=0 
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FIG. 1: The Poincare section in the phase space. 



to appear in excited-state. For simplicity, we will use the 
units of = 1 and m = 1, and we take the first excited- 
state for example. First we show the Poincare section in 
the phase space in Fig. 1 by solving the coupled equa- 
tions (9.5a) and (9.5b). 

From Fig. 1 , we can easily select the initial conditions 
as the following: 

Chaos State: A = 0.0,1 = 1.48565, G = 0.35, G = 
0.0; e = 1.0; 

Regular State: A = 0.0, A = 1.17969, G = 0.225, G = 
0.0; e = 1.0; 

And we calculated their Lyapunov Exponents (LE) to 
verify our selection. We find the LE in the first condition 
is positive and the second one is zero. Fig. 2 shows the 
FS of the probability density in the first excited state 
(n=l). 
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FIG. 2: Fourier spectra of the time evolution of the probabil- 
ity, \(j>n {x, G {t))f, where n=l. 

Prom Fig. 2., we can see that in the first excited-state, 
those 'fundamental frequencies ' still can be found in the 
Fourier spectra, such as i/f^=0.4356 Hz, i/f =0.3304 Hz, 
etc. They still play important roles in combining the 
frequencies in the FS. 

But if we compare the FS in the first-excited state and 



the ground states, we can find much difference between 
them. 

(a) . It's obvious that in the first-excited state the FS 

in regular state appears more regular than the one in 
ground state, while the FS of chaotic state appears more 
chaotic than the one in the ground state. 

(b) . We can find that in the first-excited state the 
fundamental frequencies in regular state are larger than 
the corresponding ones in the ground state, while the 
fundamental frequencies in chaotic state are smaller than 
the corresponding ones in the ground state. 

In fact, these phenomena can be explained as follows: 

(a) . It is pointed ([8]) that the 'fundamental frequen- 
cies' originate only from the coupling between the clas- 
sical and quantum parts. Because of that reason, the 
'fundamental frequencies ' will also appear in the excited 
states. 

(b) . In the first-excited state, the LE of chaotic state 
(G (0) = 0.35) is almost two times larger than that in the 
ground state (G (0) = 0.5). And we find the LE in regular 
state (G (0) = 0.225) converges to zero faster than that 
in the ground state (G(0) = 0.35). So in excited-state, 
the chaotic state will become more chaotic, and regular 
state will appear more regular. 

From above discussion, we can get a conclusion that 
the difference between regular and chaotic state in FS 
still clearly exists in the first excited state. Because of 
that reason we can use the FS to distinguish the chaotic 
and regular states even in excited-states, too. 

Now, we study the correlation function in the excited 
states. The Correlation Function (CF) is defined as 

+L +00 

CF{t)= lim ^ [ dt' [ dx\(l){x,t + t')\\(l){x,t')\, 

L—>oo 2L J J 
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(17) 

We have got the excited wave function (Eq. 11), so we 
can get the time-correlation functions in the chaotic and 
regular state and their FS, respectively. 
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FIG. 3: Fourier spectra of the time-correlation functions of 
\4>n {x,G , n=l. a) in chaotic state; b) in regular state. 



FIG. 4: Sensitivity of the probability density to the initial 
conditions for different energies, a) Ei = 2.7; b) Ei = 3.0. 
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In Fig. 3, we give out the Fourier spectra of time- 
^ correlation function in the first excited state, in which the 
^rncan value has been subtracted from the FS to remove 
tnfe big dc-component in the FS. In these figures, we can 
nnd that those 'fundamental frequencies ' appear again, 
and these figures can show the 'fundamental frequencies ' 
more clearly. So we know that, the FS of time-correlation 
can be a justifiable way to distinguish the regular and 
chaotic states. And when the system is in the chaotic 
state, the FS of the correlation function becomes desul- 
tory, and it looks to have a lot of 'noise' components. 

Finally, wc examine the sensitivity of the cigcnfunc- 
tions in the chaotic state to the initial conditions, which 
is usually considered as a basic characteristic of chaos. 

Fig. 4 shows related numerical results of the probability 
density |0„ (a;, G vs. t for the case of Ei = 2.7 

(Fig. 4a) and 3.0 (Fig. 4b). The initial conditions are 
selected by following rule: all parameters are the same 
except that G(0) has very little difference between two 
solutions, which are G(0) = 0.35 and G(0) = 0.3501. We 
can sec that in the first excited state, there still exists 
such sensitivity for the eigenfunction, which is influenced 
by the total system energy Ei. And according to Ref. 
[8], when the system energy rises, the chaotic behavior 
becomes stronger. 
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IV. CONCLUSIONS 

In this paper, we analytically calculated the excited- 
state wave functions of a quantum oscillator coupled with 
a classical harmonic oscillator. Instead of using the def- 
inition of the classical chaos, such as the sensitivity of 
the system to the initial conditions, we have found the 



'fundamental frequencies ' in the ground state still appear 
in the excited state, and found a very useful way to get 
these characteristics. We think these characteristics can 
be used to distinguish the regular and chaotic states when 
the so-called semiquantum chaos emerges in the actual 
system and our result is helpful to further research work 
in this field. 
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